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ABSTRACT 

The shape of the probability distribution function (PDF) of molecular clouds is an important in¬ 
gredient for modern theories of star formation and turbulence. Recently, several studies have pointed 
out observational difficulties with constraining the low column density (i.e. Ay < 1) PDF using dust 
tracers. In order to constrain the shape and properties of the low column density probability distribu¬ 
tion function, we investigate the PDF of multiphase atomic gas in the Perseus molecular cloud using 
opacity-corrected GALFA-HI data and compare the PDF shape and properties to the total gas PDF 
and the N(H 2 ) PDF. We find that the shape of the PDF in the atomic medium of Perseus is well 
described by a lognormal distribution, and not by a power-law or bimodal distribution. The peak of 
the atomic gas PDF in and around Perseus lies at the HI-H 2 transition column density for this cloud, 
past which the N(H 2 ) PDF takes on a powerlaw form. We find that the PDF of the atomic gas is 
narrow and at column densities larger than the HI-H 2 transition the HI rapidly depletes, suggesting 
that the HI PDF may be used to find the HI-H 2 transition column density. We also calculate the 
sonic Mach number of the atomic gas by using HI absorption line data, which yields a median value 
of Ms = 4.0 for the CNM, while the HI emission PDF, which traces both the WNM and CNM, has a 
width more consistent with transonic turbulence. 

1. INTRODUCTION 

Understanding the process of star formation is compli¬ 
cated by the interaction of gravity, magnetic fields, stellar 
feedback, and turbulence. The details of the collapse of 
molecular clouds on parsec scales may determine the key 
properties of the initial mass function and the star for¬ 
mation rates in galaxies (e.g. see Hennebelle & Falgarone 
2012). Thus in order to obtain a complete picture of the 
process of star formation and the evolution of galaxies, 
the dynamics of molecular clouds must be understood. 

Towards this end, the column density probability dis¬ 
tribution function (PDF) is a common tool used by theo¬ 
rists and observers alike to disentangle the role of turbu¬ 
lent motions, stellar feedback, gravitational contraction 
and magnetic support. Several models claim that the 
high density end of the PDF should provide the dense 
gas fraction (Krumholz & Mckee 2005; Hennebelle & 

Chabrier 2011; Padoan & Nordlund 2011). Furthermore, 
the highest column density regime of the PDF of self- 
gravitating turbulent clouds has a power-law distribution 
as demonstrated in numerical simulations (Federrath & 

Klessen 2012, 2013; Collins et al. 2012; Burkhart, Collins 
& Lazarian 2015) and observations (Kainulainen et al. 

2009; Froebrich & Rowles 2010; Schneider et al. 2013, 

2105; Lombardi, Alves & Lada 2015). 

The lower column density material in the PDF 
of molecular clouds should trace primarily non-self- 
gravitating gas, be dominated by turbulent motions, 
and is expected to have a lognormal form (Vazquez- 
Semadeni 1994; Padoan et al. 1997; Burkhart & Lazar¬ 
ian 2012). Lognormal PDFs have been observed in mul- 
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tiple phases in the ISM including CO (Brunt 2010), in 
dust (e.g., Kainulainen & Tan 2013; Froebrich & Rowles 
2010; Schneider et al. 2013; Alves de Oliveira et al. 
2014), and in the diffuse ISM (Berkhuijsen & Fletcher 
2008; Burkhart et al. 2010; Burkhart, Gaensler, Lazar¬ 
ian 2012; lacobelli et al. 2014). The PDF has also been 
noted by several authors to be an an important method 
for determining the sonic Mach number {Ms, Federrath 
et al. 2008; Burkhart et al. 2009; Kainulainen & Tan 
2013; Burkhart & Lazarian 2012). However Burkhart, 
Collins, & Lazarian (2015) have reported that the PDF- 
sonic Mach number relation should only be robust for 
diffuse gas and not in a self-gravitating medium. 

Tracing the PDF from observations is non-trivial as 
multiple tracers must be used in order to sample the 
largest dynamic range of densities. For dense material 
in molecular clouds, dust emission and absorption are 
the most promising tracers of the true column density 
distribution with the largest dynamic range of densities. 
This is in contrast to molecular line tracers such as CO, 
which suffer from opacity, excitation temperature, ab¬ 
sorption and depletion effects (Goodman et al. 2009; 
Burkhart et al. 2013a,b). However, recently several au¬ 
thors (e.g. Schneider et al. 2015; Lombardi, Alves, & 
Lada 2015) have noted that dust emission and extinction 
are problematic probes of the low column density mate¬ 
rial in molecular clouds. The PDF as measured from 
observations of dust can suffer several biases including 
resolution, noise, boundary effects and line-of-sight con¬ 
tamination from the foreground and background. Partic¬ 
ularly, Lombardi, Alves & Lada (2015) pointed out that 
the lognormal portion of the PDF (i.e. the portion of the 
PDF traced below A^ < 0.1 mag or Ay < 1 mag) may 
not be securely traced by dust. In light of this, other 
tracers should be employed to determine the properties 
of the PDF in molecular clouds at lower densities. 

Our motivation for this Letter is to explore the use of 
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HI as a tracer of the low column density PDF in star 
forming molecular clouds. In general, HI is an excel¬ 
lent tracer of low column density atomic gas in the far 
outskirts of molecular clouds. While it can be optically 
thick closer to GMC centers (see Bihr et al. 2015), HI 
absorption measurements can be obtained to provide full 
column density information and provides insight into the 
nature of the atomic to molecular transition in GMCs 
(Lee et al. 2015). In particular we focus our study on 
the Perseus molecular cloud, whose atomic and molecu¬ 
lar gas and dust content, as well as the HI-H 2 transition, 
has been well-studied in a number of previous works (Lee 
et al. 2012; Lee et al. 2014; Stanimirovic et al. 2014; 
Lee et al. 2015). 

This Letter is organized as follows: in Section 2 we 
provide a description of the HI and dust data sets of the 
Perseus molecular cloud used in this study. In Section 3 
we discuss our findings regarding the PDF of HI in and 
around Perseus and compare it to the dust PDF. In Sec¬ 
tion 4 we discuss our results followed by our conclusions 
in Section 5. 

2. DATA 

In this section, we provide a summary of the data in 
our study and refer to Lee et al. (2012; 2015) and Stan¬ 
imirovic et al. (2014) for more details. 

2.1. HI Emission and Absorption Data 

We use the HI data cube from the GALFA-HI Sur¬ 
vey (Peek et al. 2011), centered at (R.A.,decl.) = 
(03‘'29“52",-b30°34'l") in J2000 with a size of ~15° x 
9°. To derive the HI column density, Lee et al. (2012) 
integrated the HI emission from Vlsr = — 5 km s“^ to 
-|-15 km s“^ under the optically thin assumption. This 
HI velocity range for Perseus was determined based on 
the maximum correlation between the derived A^(HI) and 
the 2MASS Ay from the COMPLETE Survey (Ridge et 
al. 2006). To construct an accurate PDF of atomic gas, 
we use the opacity-corrected HI column density image 
(see Lee et al. 2015 for details), but exclude the regions 
with A^(HI) < 6x 10^° cm“^ from our analysis due to 
possible contamination with the Taurus molecular cloud 
(see Figure 1). 

We also use HI absorption data from Stanimirovic et 
al. (2014) in this work to estimate the sonic Mach num¬ 
ber in the CNM of Perseus. The HI emission-absorption 
spectral pairs were obtained toward 26 radio continuum 
sources located behind Perseus, and the Gaussian de¬ 
composition method (Heiles & Troland 2003a) was em¬ 
ployed to calculate the properties of individual CNM and 
WNM components (e.g. spin temperature, optical depth, 
“true” HI column density, etc.). 

2.2. Ay based on Far-infrared Dust Emission Data 

We use the total gas (hereafter referred to as the Ay 
map) and N(H 2 ) images from Lee et al. (2015) and refer 
to that work for details on the data processing. The 
N(H 2 ) image is created from the iV(HI) and Ay images 
as: 

N{H,) = ^i^-N{HI)) ( 1 ) 

using the local dust-to-gas ratio (D/G) of ~1 x 10“^^ 
mag cm^ from Lee et al. (2015). This is D/G value we 


adopt for the comparison between the HI and dust PDFs 
(Section 3). In this case, the N(H 2 ) PDF represents the 
molecular gas without the atomic contribution. The final 
Ay image (Figure 1) is at 4.3' angular resolution with a 
pixel size of 4.3' (same for the N(HI) image), and the 
median Icr uncertainty in Ay is ~0.17 mag. 

3. RESULTS 

We show the PDF of atomic and molecular gas associ¬ 
ated with the Perseus molecular cloud in Figure 2. The 
black line shows the PDF of the Ay image from Lee et al. 
(2015), which peaks close to Ay = 1 mag (straight verti¬ 
cal line), which is the column density threshold Lee-|-15 
found for the HI-to-H2 transition in Perseus (tan box), 
as well as the Ay limit Lombari, Alves & Lada (2015) 
suggested for dust as a unreliable tracer of the low col¬ 
umn density PDF. The red dashed line shows the PDF 
of atomic gas. The blue line shows the N(H 2 ) PDF from 
Lee et al. (2015). The PDFs as shown in Figure 2 were 
binned into 100 bins and have been normalized by the to¬ 
tal number of points in the images (which are the same 
in both the HI image and the dust images). The skew¬ 
ness and kurtosis of the HI GALFA data is log(HI), 0.1 
and -0.05 respectively, which suggest the distribution is 
nearly Gaussian and is within the range expected from 
MHD simulations with no gravity (i.e. see Burkhart, 
Collins & Lazarian 2015, their Figure 6). We fit a Gaus¬ 
sian distribution to the HI PDF and the Ay (black lines 
downshifted for ease of comparison with the shape) and 
Hnd the Htted standard deviation to be a hi = 0.13 and 
<^Av = 0.11. The 1-sigma errors in the HI image are 0.08 
mag (when D/G = 1 x 10“^^ is adopted). We overplot 
Poisson error bars for both the Ay and HI PDFs. 

The HI gas associated with Perseus can be well de¬ 
scribed by a narrow lognormal distribution which peaks 
close to Ay = 1 mag, which we mark with a straight 
vertical line. The shape of the N(H 2 ) PDF at column 
densities above the HI-H 2 threshold is dominated by a 
power law tail, similar to the Ay PDF, while at low col¬ 
umn densities below the HI-H 2 threshold the N(H 2 ) PDF 
flattens due to the fact that there is a large amount of 
Hl-dominated diffuse gas around Perseus, which is sub¬ 
tracted from the low amount of total gas via Equation 1. 
The shape of the total gas Ay PDE (which includes both 
N(H 2 ) and HI gas) below the HI-H 2 threshold is domi¬ 
nated by a lognormal shape (i.e. the HI PDF) while 
above it displays a powerlaw tail. The fact that the 
N(H 2 ) pdf is not lognormal at column densities below 
the HI-H 2 threshold suggests that the lognormal portion 
of the Ay PDF is produced by the HI gas, which of 
course is removed to create the N(H 2 ) PDF. The lognor¬ 
mal behavior of the atomic gas is in agreement with nu¬ 
merical simulations of supersonic MHD turbulence (e.g. 
Burkhart & Lazarian 2012; Federrath & Klessen 2012; 
Burkhart Gollins & Lazarian 2015) and only the dense 
self-gravitating molecular gas displays a power-law tail 
behavior. 

We note that the PDF of the HI column density distri¬ 
bution for Perseus is narrow, which is expected consider¬ 
ing the very uniform distribution of HI in this molecular 
cloud (see Lee et al. 2012; 2015). This is because hy¬ 
drogen is converted from atomic to molecular (HI-to-H 2 
transition) at A^(HI) ~ (9“11) x 10^° cm“^ for Perseus 
(Lee et al. 2015). However, in general the uniform HI 
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Fig. 1. — Top: Opacity-corrected HI column density image from Lee et al. (2015). Bottom: image from Lee et al. (2015). The blank 

pixels correspond to the regions with possible contamination, e.g. point sources and a background HII region. In addition, the regions with 
N(HI) < 6 X 10^® cm~^ are shown with black crosses and are not used for our analysis as these pixels could be associated with the Taurus 
molecular cloud. See Lee et al. (2012) for details on the pixel masking. 
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Fig. 2.— The PDF of atomic and molecular gas in the Perseus cloud using the Av and N(H2) images from Lee et al. 2015 (black and 
blue dashed lines, respectively) and the GALFA 21cm emission (red dashed line), as discussed in the text. The dotted black line directly 
under the red line is a Gaussian fit to the HI PDF, downshifted for ease of visibility. The PDFs are normalized by the total number of 
data points, which are identical over both data sets. The tan box brackets the range of densities for the HI-H 2 transition as derived in Lee 
et al. (2015). The solid black vertical line is a references for Ay = !)• The median l-a error level for the HI data is 0.08 mag (when D/G 
= 1 X 10“^^ is adopted) and for the Ay data 0.17 mag. 
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column density distribution could also result from a large 
amount of optically thick HI gas (see Bihr et al. 2015). 
The amount of the optically thick HI in Perseus was stud¬ 
ied in Lee et al. (2015). They found that the contribution 
from the optically thick HI is minor (less than 20%) and 
H 2 formation is indeed responsible for the nearly uni¬ 
form HI column density across Perseus (see also Bialy 
et al. 2015). The fact that the HI distribution is not 
bimoda0 suggests that much of the gas may be close the 
to transition to molecular gas. 

Interestingly, the HI and Ay PDFs peak at comparable 
column density values (within the 1-a errors of the data) 
and have similar values for the fitted Gaussian widths, 
however a number of effects must also be considered. For 
the HI PDF, the peak value will depend on the D/G 
that is used as well as the cloud superpositions along the 
line of sight, as well as the amount of cold HI. These 
effects have been well-studied in Perseus and the line- 
of-sight velocity which is used to disentangle different 
clouds, as well as correction for high optical depth, is well 
constrained (see Lee et al. 2015). For the dust PDF, the 
peak location depends on the LOS contamination, offset 
corrections and superposition effects (see Schneider et 
al. 2015; Lombardi, Alves & Lada 2015). Stanchev et 
al. (2015) demonstrated different selected regions in the 
Perseus PDF of Planck dust optical depth at 353 GHz 
(^■ 353 ) would show different values of the PDF peak and 
width. This is most likely the reason why the peak of 
the Ay PDF for Perseus in the Lee et al. (2015) data 
is offset from the peaks of the Kainulainen et al. (2009) 
and Sadavoy et al. (2014) Perseus PDFs (Ay 0.8 
as compared to Ay ss 1.5, respectively). The width of 
the Perseus PDF reported in Kainulainen et al. (2009) 
is much larger (cr = 0.47) than the width found in our 
work, which is more compatible to the value reported by 
Stanchev et al. (2015) using Planck data, i.e. a « 0.11. 

3.1. The Sonic Mach Number 

Numerical simulations of MHD turbulence have sug¬ 
gested that the sonic Mach number is related to the width 
of the lognormal PDF (e.g. see Federrath et al. 2008, 
Burkhart & Lazarian 2012 and references therein). We 
compute the sonic Mach of the HI gas of Perseus using 
both the PDF and GNM absorption measurements. 

First, we compute the sonic Mach number using the 
absorption line data of Stanimirovic et al. (2014) and 
Equation 17 of Heiles & Troland (2003): 

m 2 = 4.2(Zi^ - 1) (2) 

We show a histogram of the Mach numbers derived 
from Equation 1 using 107 GNM absorption components 
(described in Section 2 and in Stanimirovic et al. 2014) in 
Eigure 3. The median value is Mg =4.0, and is indicated 
by the vertical dashed line. 

Since the HI distribution shows a single lognormal 
distribution we can estimate the sonic Mach number 
from the width of the lognormal using Equation 4 from 
Burkhart & Lazarian (2012): 

Bimodal distributions of HI are also predicted due to the fact 
that this tracer probes a multiphase medium with regimes of ther¬ 
mal instability, see Gazol & Kim (2013). 



Fig. 3.— Histogram of the sonic Mach number as calculated from 
the absorption line data for Perseus. The median value is Ms = 4.0 
and is shown with a straight vertical line. 


c^Le/Eo = ^ ^ Inih^Ml + 1) (3) 

for different values of parameters A (which depends 
on the optical depth, see Burkhart et al. 2013a) and b 
(which depends on the type of turbulence driving, see 
Federrath et al. 2008). Figure 4 shows curves of Equa¬ 
tion 2 for different values of A and b. The narrow PDF 
of the HI gas in Perseus is consistent with transsonic 
Mach numbers or with very low values of the parameter 
A. Low values of A were shown to be consistent with 
a more optically thick medium (Burkhart et al. 2013a), 
therefore it is more likely that the PDF of the HI gas, 
which is tracing a mix of the WNM and GNM, has a 
lower average sonic Mach number than the sonic Mach 
number reported from absorption lines, which only trace 
the GNM. 
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Fig. 4.— Variance- sonic Mach number relationship given by 
Equation 2 for different values of A and b. We use the same distri¬ 
bution of sonic Mach number as presented in Figure 3. 
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4. DISCUSSION 

The shape and properties of the PDF of molecular 
clouds are important ingredients for modern theories of 
star formation. In this letter we found that the low col¬ 
umn density material in the Perseus molecular cloud as 
traced by neutral hydrogen takes on a lognormal shape 
rather than a bimodal or power-law distribution. A num¬ 
ber of processes can induce a lognormal distribution in 
density, including subsonic and supersonic turbulence 
but also gravitationally contracting clouds which are not 
necessarily turbulent (see Tassis et al. 2010). However, 
high Reynolds number flows (with Re = LVjv > 10^) 
are almost always turbulen10 and the density of HI is 
diffuse enough to ignore the effects of self-gravity. If tur¬ 
bulence is indeed responsible for the lognormal HI PDF, 
our study along with Lee et al. (2012, 2015) and Bialy 
et al. (2015) suggests that the chemical and dynam¬ 
ical evolution of the interstellar medium can be mod¬ 
eled separately, at least for low-mass molecular clouds 
like Perseus. In Lee et al. (2012, 2015) and Bialy et 
al. (2015), steady state and chemical equilibrium mod¬ 
els of H 2 formation (Krumholz et al. 2009; Sternberg et 
al. 2014) were tested against the observations of Perseus 
and a good agreement was found. This good agreement 
implies that, at least for clouds where strong stellar feed¬ 
back is not important, H 2 formation can be modeled in 
steady state while turbulence determines the density dis¬ 
tribution of the cloud. 

In this letter we derived the sonic Mach number of the 
atomic gas in and around Perseus using two independent 
methods: the width of the PDF of the HI column density 
and absorption line fitting using data from Stanimirovic 
et al. (2014). The median Mach number derived from the 
absorption line data {Mg = 4.0). On the other hand, the 
width of the HI PDF indicates that HI (CNM -|- WNM) 
in and around Perseus is on average transsonic with the 
sonic Mach number lower than what the HI absorption 
line data and CO observations suggest (e.g. Padoan et al. 
2003; Ms > 10 for Perseus based on CO observations). 
The 21-SPONGE survey (Begum et al. 2010; Murray et 
al. 2014, 2015) will provide additional absorption line 
data using the VLA and will be able to further test the 
HI PDF-Mach number relation for other clouds. Ulti¬ 
mately knowing the Mach number in the atomic gas will 
constrain the initial conditions for numerical simulations 
of diffuse and star forming clouds. 

It should be pointed out that there are observational 
difficulties to derive accurate HI column densities asso¬ 
ciated with molecular clouds. In particular one must 
investigate clouds that lie off the plane of the Galaxy 
to mitigate line-of-sight confusion (e.g. Perseus at the 
Galactic latitude of —20°). It is also important to define 
the HI velocity range for an individual molecular cloud, 
correct for high optical depth HI, and check that the HI 
correlates to some degree with the dust map. All of these 
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procedures were performed for the data set used in this 
study (see Lee et al. 2012; 2015). 

Despite the difficulties in accurately constraining the 
PDF of molecular clouds through HI, the atomic gas PDF 
should be compared to the molecular PDFs (i.e. traced 
by dust emission/absorption) whenever possible as they 
probe different column density regimes. The peak of both 
the Av and HI PDFs explored in this paper lies roughly 
at the HI-to-H 2 transitional column density found in Lee 
et al. (2012) for Perseus. The HI PDF peak in GMCs 
may trace this transitional density since the HI rapidly 
depletes at higher column densities to H 2 . Future studies 
should investigate the HI column density distribution in 
other GMCs and compare with the measured HI-to-H 2 
transition as our study suggests the transition between 
the lognormal to power-law behavior of the PDF traces 
the HI-to-H 2 transition. 

5. CONCLUSIONS 

We investigate the multiphase PDF of the atomic gas 
in the Perseus molecular cloud and compare it to the 
PDF traced by dust. We find that: 

• The shape of the PDF of the atomic gas in Perseus, 
which traces gas with Ay < 1 mag, is well de¬ 
scribed by a lognormal distribution, and not by a 
power-law distribution or a bimodal distribution. 

• The shape of the N(H 2 ) PDF is a power law past 
the HI-to-H 2 transition, similar to the total gas 
(Av) PDF. 

• The transition from the lognormal atomic PDF to 
the power-law dust PDF happens roughly at the 
HI-to-H 2 transition. This suggests that the HI and 
dust PDFs may be used in other clouds to find the 
HI-to-H 2 transition. 

• We calculate the sonic Mach number of the atomic 
gas by using HI absorption line data, which yields 
a median value of Mg = 4.0 for the CNM while the 
HI emission PDF, which traces both the WNM and 
CNM, has a width more consistent with transonic 
turbulence. 
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